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Abstract We present solutions of the Dirac equation with spin symmetry for vector and scalar 
modified Poschl- Teller potential within framework of an approximation of the centrifugal term. 
The relativistic energy spectrum is obtained using the Nikiforov-Uvarov method and the two- 
component spinor wavefunctions are obtain are in terms of the Jacobi polynomials. It is found 
that there exist only positive-energy states for bound states under spin symmetry, and the 
energy levels increase with the dimension and the potential range parameter a. 
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1- 

Solutions to relativistic equations play a very important in many aspects of mordern 
physics. In particular, the Dirac equation which describe the motion of a spin-^ particle has 
been used in solving many problems of nuclear aand high-energy physics. The spin symmetry 
arises if the magnitude of the attractive Lorentz scalar potential S(r) and the replusive vector 
potential V(r) are nearly equal, i.e S(r) ~ V(r), in the nuclei, while the pseudo-spin symmetry 
occur when S(r) ~ —V(r) [1-3]. The case of the exact spin and pseudo-spin symmetries has 
been shown to correspond to the SU(2) symmetries of the Dirac Hamiltonian [3]. The spin 
symmetry is relevant for mesons [4] and the pseudo-spin symmetry is used to explain deformed 
nuclei [5], super-deformation [6] and to establish an effective nuclear shell-model scheme [7- 
9]. Also, various potentials such as the Morse potential [10-12], Wood-saxon potential [13], 
Coulomb and Hartmann potentials [14], Eckart potential [15, 16], Poschl- Teller potential [17, 
18] and the harmonic potential [19, 20] have been studied within the frame work of the spin 
and pseudospin symmetries. 
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Moreover, with the interest in higher dimensional theory, the multi-dimensional quantum 
mechanical equations-relativistic and nonrelativistic- have been solved with various physi- 
cal potentials. To mention a few, the D-dimensional Schrodinger has been studied with the 
Coulomb- like potential [21], pseudoharmonic potential [22], Hulthen potential [23] and Poschl- 
Teller potential [24]. In addition, the D-dimensional relativistic Klein-Gordon and Dirac equa- 
tions have been studied with many exactly solvable models [25-30]. However, some physical 
potentials can only be solve exactly for the s-states. Unfortunately, the modified Poschl- Teller 
potential is one of these potentials. For instance, a recent work [31] has presented the s-wave 
solutions of the Dirac equation with the Poschl- Teller potential under the conditions of the 
exact spin symmetry and pseudospin symmetry. However, in order to extend the solutions of 
the modified Poschl- Teller potential to any £ ^ state, some recent studies [17, 24] have used 
a hyperbolic approximation for the centrifugal term to obtain the non-relativistic solutions 
of the modified Poschl- Teller potential. In light of this, the present Letter intends to extend 
the discussions on the relativistic Poschl- Teller potential to D-dimensions by presenting the 
bound-state solutions of the D-dimensions Dirac equation with spin symmetry for Lorentz 
vector and scalar modified Poschl- Teller potential using the Nikiforov-Uvarov method [32]. 

The D-dimensional Dirac equation with a scalar potential V s (r) and a vector potential 
V v (r) and mass /j, can be written in natural units H = c = 1 as [27, 33, 34] 

D 

H*(r) = E nrK ^(r) where H = £ a jPj + /% + V s (r)} + V v (r) (1) 

3=1 

where E„ rK is the relativistic energy, {<x/} and $ are Dirac matrices, which satisfy anti- 
commutation relations 

ajdik + akdij = 25jkl 

atjfl + paj = (2) 
a 2 = ? 2 ' =1 

and 

d 

Pj = ~ id j = ~ i 'Q^T. 1 ^ D - ( 3 ) 

The orbital angular momentum operators Ljk, the spinor opertaors Sjk and the total angular 
momentum operators Jjk can be defined as follows: 

d . d 

Lj k = —L jk = ixj-g— — ixk-g^-, Sjk = —Skj = i&j6tk/2, Jjk = Lj k + Sjk- 

D D D 

L 2 = Y,L%, S 2 = Y,S] k , J 2 = Y, J h> l^j<k^D. (4) 

j<k j<k j<k 

For a spherically symmetric potential, total angular momentum operator Jj k and the spin- 
orbit operator K = —(3(J 2 — L 2 — S 2 + (D — l)/2) commutate with the Dirac Hamiltonian. 
For a given total angular momentum j, the eigenvalues of K are k = ±(j + (D — 2)/2); 
k = — (j + (D — 2)/2) for aligned spin j = £ + \ and k = (j + (D — 2)/2) for unaligned spin 

Thus, we can introduce the hyper spherical coordinates [35] 
x\ = r cos 9\ 

x a = r sm6\ . . . sin# a _i cos0, 2 ^ a ^ D — 1 (5) 
xd = r sin#i . . . sin#D_2 sin0, 
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where the volume element of the configuration space is given as 

D D-l 

Y\ d Xj = r D ~ l drdtt dtt = Yl ( f anO j ) j - 1 d6 j (6) 

3=1 3=1 

with ^ r < oo, ^ 9^ ^ tt, k = 1, 2, . . . D — 2, 0^0^ 27r, such that the spinor 
wavefunctions can be classified according to the hyperradial quantum number n r and the 
spin-orbit quantum number k and can be written using the Pauli-Dirac representation 



D-l 
2 



/ F nrK (r)Yf m (n D ) \ 



V iGn rK (r)Yf m (n D ) J 



(7) 



where F nrK (r) and G nrK {r) are the radial wave function of the upper- and the lower-spinor 
components respectively, Yj m (Qr)) and Yj m (Qrj) are the hyperspherical harmonic functions 
coupled with the total angular momentum j. The orbital and the pseudo-orbital angular 
momentum quantum numbers for spin symmetry I and and pseudospin symmetry t refer to 
the upper- and lower-component respectively. 

Substituting Eq. (7) into Eq. (1), and seperating the variables we obtain the following 
coupled radial Dirac equation for the spinor components: 



^ + 1 ) F "<-"- ( . r) ~ [/' + /; <-' -Mr)]O lh ,Ar) 



d_ 

dr 



GnUr) = \m~ E n T n + S(r)]F„ rK (r) 



(8) 



(9) 



where A(r) = V v (r) - V s (r), S(r) = V v (r) + V s (r) and k = ±(2£ + D- l)/2. Further details 
of the derivation can be obtain from refs [36-38]. Using Eq. (8) as the upper component and 
substituting into Eq. (9), we obtain the follwoing second order differential equations 



d 2 k(k + 1) 
dr 2 r 2 



[JJ. + E nrK - A(r)][/z - E nrK + E(r)] + 



dA(r) /_d_ 
dr V dr 



dr 2 



k(k — 1) 



- [fi + E nrK - A(r)][/i - E nrK + E(r)] 



Mr) + E nrK - A(r)] 



dr 



- E nrK + S(r) 



Fn rK (r) = 



(10) 

G nrK {r) = 



(11) 

We note that the energy eigenvalues in these equation depend on the angular momentum 
quantum number I and dimension D. However, to solve these equations, we shall use an 
approximation for the centrifugal barrier and obtain the solutions using the Nikiforov-Uvarov 
method. 

Next, we give a brief description of the conventional Nikiforov-Uvarov method. A more 
detailed description of the method can be obtained the following reference [32]. With an 
appropriate transformation s = s(r), the one dimensional Schrodinger equation can be reduced 
to a generalized equation of hyper geometric type which can be written as follows: 



a(s) 



a 2 (s) 







(12) 



Where a(s) and a(s) are polynomials, at most second-degree, and f(s)is at most a first-order 
polynomial. To find particular solution of Eq. (12) by separation of variables, if one deals with 



ip(s) = (f)(s)y nr (s), 



(13) 
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Eq. (12) becomes 

a(s)y'l+r{s)y' nr + \y nr =Q (14) 

where 

*(.) = »(•>$$ as) 

r( s ) = f(s) + 2it(s),t'(s) < 0, (16) 



7r(s) 



u — r 



a — t 



a + ta, 



and 



A = t + ir'(s). 



(17) 



(18) 



The polynomial r(s) with the parameter s and prime factors show the differentials at first 
degree be negative. However, determination of parameter t is the essential point in the 
calculation of ir(s). It is simply defined by setting the discriminate of the square root to zero 
[32]. Therefore, one gets a general quadratic equation for t. The values of t can be used for 
calculation of energy eigenvalues using the following equation 



\ j. i i( \ U \ n r\ n r 1) II / \ 

A = t + 7r (s) = —n r T (s) a (s). 



(19) 



Furthermore, the other part y nr (s) of the wave function in Eq. (12) is the hypergeometric-type 
function whose polynomial solutions are given by Rodrigues relation: 



Vn r {s) 



J->n r 11 

p(s) ds nr 



W nr (s)p(s)} 



(20) 



where B rir is a normalizing constant and the weight function p(s) must satisfy the condition 
[32] 

W = rp. (21) 

The Lorentz vector V v (r) and scalar V s (r) modified Poschl- Teller potential can be defined 
as follows [24, 39-41] 



V v (r) 



cosh (or) 



and V s (r) 



So 



cosh (ar) 



(22) 



where a is related to the range of the potential and Vq and So are the depths of the vector and 
scalar potentials respectively. Moreover, we can approximate the centrifugal terms as follows 
[17, 24] 



1 



a 



sinh (ar) 



Substituting Eqs. (22) and (23) into Eqs. (10) and (11), we have 



d 2 ci 2 k(k + 1) 



dr 2 sinh 2 (ar) 

d 2 o?k{k — 1) 
dr 2 sinh 2 (ar) 



[p + E nrK - A(r)][p - E nrK + E(r)] + 



dA(r) 
dr 



(dr + r) 



Mr) + E nrK - A(r)] 



[p + E nrK - A(r)][p - E nrK + E(r)] 



dr \dr r ) 

W) - E nrK + E(r) 



(23) 



Fn r n{r) = 

(24) 
G nrK (r) = 
(25) 
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where 



A(r) 



So-V 



and S(r) 



(26) 



cosh 2 (or) cosh 2 (or) 

For the case of spin symmetry, V v (r) ~ V^(r), i.e. A(r) = K,(r) — V,(r) = C\ (a constant), 
which implies that dA }^ = 0. Thus, putting this into Eq. (24), we have 



d 2 a 2 K( K + l) (y + s )(£ nrK + /U _c 1 ) 
(/i - E nrK )\^l + £/„ rK - UiJ H 



_dr 2 sinh 2 (ar) 



cosh (or) 



If we take the transformation s = tanh 2 (ar), Eq. (27) becomes 



F nrK (r) = 0. 

(27) 



4s 2 (l 



:[-<5 S 2 + (<5 + 7-e 2 ) S - 7 ]F nr(t = (28) 



where 



6 2 = 



(//- E nrK )(n + E nrK - d) 



5 = 



(V + S )(E nrK + »-d) 



and 7 = k(k + 1). (29) 



or or 
Comparing Eqs. (28) and (12) we can define the following 

f(s) = 1 - 3s, a(s) = 2s(l - s) and a(s) = -5s 2 + (7 + 6 - e 2 )s - 7 (30) 

Inserting these into Eq. (17), we have the following function 



ir(s) 



1 



^ ± i 7(1 + 45 " 8 *> 2 + ( 8t - 4(7 + ^ - e 2 ) - 2)s + 4 7 + 1 (31) 



The constant parameter i can be found by the condition that the discriminant of the expression 
under the square root has a double root, i.e., its discriminant is zero. Thus the possible value 
function for each value of t is given as 



tt(s) 



1 



\ [(-2e + 7TT47~)s-7TT47~] for t = -±(7 - 5 + e 2 ) + \e^T+A^ 

\ [(2e + 7T+47>- VTT4t~] for t = -±(7 - 8 + e 2 ) - \e^T+A^ 

(32) 

By Nikiforov-Uvarov method, we made an appropriate choice of the function 7r(s) = — 
\ [(2e + VI + 47) s - V 1 + 47] such that by Eq. (19), we can obtain the eigenvalue equation 
to be 



1 



1 



1 



1 



- - (7 - S + e 2 ) - -e Vl+47 " 2 ( 2e + V 1 + 47) - - = n r [4 + 2e + y/l + 4 7 ] + 2n r (n r - 1) (33) 



Eq. (33) can be written in the powers of e as follows 



e 2 + e 



2(2n r + 1) + v 7 ! + 47j + (7 - <*) + [(1 + 2n r ) + + 
such that we can obtain 



0, 



l r 



-2(2n r + 1) - 71 + 47 + v 7 !^ 



(34) 



(35) 



from which we can obtain a rather complicated transcendental energy equation: 

„2 

2(2n r + 1) + (2k + 1) - 

a 



a 



(p-En^ifi+En^-d) = — 



2(2n r + 1) + (2k + 1) - i Va 2 + 4(V + 5 )(K rK + \i - d) 

(36) 



i 2 
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If we define a principal quantum number n = 2n r + 1 + 1 , Eq. (36) becomes 



(M-£„)0i + £„-Ci) = 



2n + D - - + 4(V + S )(E n + fi - d) 
a 



(37) 



where we have chosen k = (2£ + D — l)/2 and n = 1, 2, 3, . . . Some numerical values of the 
energy levels E(a,n,D) for some dimensions and exited states are given in Table 1. 

We now obtain the spinor components of the wavefunction for the spin symmetry case 
using the Nikiforov-Uvarov method. By substituting tt(s) and a(s) into Eq. (15), and solving 
the first order differential equation to have 



0(s) = s( K+1 v 2 (i- s y/ 2 . 

Also using Eq. (18), the weight function p(s) can be obtained as 
Substituting Eq. (39) into the Rodrigues relation (20), we have 

rl n r r 

VnM = B nr S-^-V/\l - JL_ [ s n, + (2«-l)/2 (1 _ s) n r + 

Therefore, we can write the upper component F UrK (s) as 

F nrK (s) = C nr s^'\l - Sr /2 P ((2«-DA «) ( i _ 2s) 



(38) 



(39) 



(40) 



(41) 



where C nr is the normalization constant, and we have used the definition of the Jacobi poly- 
nomials [42], given as 



(-iy 



d" 



(l-s) a+n (l + s 



n!2 n (l - s) a (l + s) b ds n 
The lower-component can be obtain as follows using Eq. (8) 

G nrK (s) = A^P^-W _ 2s) + A 2 (s)P^ 2 ' ^(1 - 2s) 

where 

C nr as^(l - s)^ (l-s)- §a] + c^-^ 



(42) 



(43) 



Ms) = 



and 



Ms) 



with constant D nr defined by 



D„ 



At + E nrK — C\ 

Ay«S (K+2)/2 (l-s) (£+2)/2 

2e + 2n r + 2k + 1 
; X ^n„ 



(44) 



(45) 



. Moreover, to compute the normalization constant C nr , it is easy to show that 



/ r^s— F nr «(r) r D - 1 dr= / |F„„ K (r)| 2 dr = / \F nrK (s)\< 
Jo Jo Jo 



ds 



2ay/s{\ - s) 



1 (46) 



where we have also used the substitution s = tanh 2 (ar). Putting Eq. (41) into Eq. (46) and 
using the following definition of the Jacobi polynomial [42] 



Pt b) («) = ^IT^ 2i?1 (" n ' a + b + n + 1; 1 + 



1 - s 



(47) 
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we arrived at 



Cn r N nr f s K+1 2(l - S y- l [ 2 F 1 (-n r ,K + e + n r + 1/2; k + 1/2; s)} 2 ds = a (48) 
J o 



where N nr = ± 



r(n r +K+l/2) 2 
n r !r(ft+l/2) 



and 2F1 is the hypergeometric function. Using the following 



series representation of the hypergeometric function 

o 

(ci)„...(c 9 ) n n\ 



rp f . \ _ \^ (oi)n-(Op)n S n , . 

p-fgiai) ■•■) %>) Cl, ...,C 9 , Sj — 2_^l r \ < r \ „! V 4y J 



ra=0 

we have 



" ^ («+l/2)ii! (- + l/2),j! Jo ( S) " " 

(50) 

Hence, by the definition of the Beta function, Eq. (43) becomes 

Using the relations B(x,y) = an( ^ *he Pocbhammer symbol (a) n = F p^"' ) , Eq. (51) 

can be written as 

r i at (-"r)fc(« + e + n r + l/2)i(« + §), ^ (-n r )j(«; + e + n r + l/2)j(n + i + §)j a 



4=0 + « + §)*(« + l/2)ifc! (e + « + *+!),•(« + 1/2)^1 S(«+f,e) 

(52) 

Lastly, Eq. (52) can be used to compute the normalization constants for n r = 0,1,2,... In 
particular for the ground state, i.e n r = 0, we have 



Co = 4/ 777 — 7- ( 53 ) 

In conclusion, the solutions of the Dirac equation with spin symmetry for the modified 
Poschl- Teller potential has been extended to a multi-dimensional case. The energy levels and 
the spinor-components of the wavefunction were obtained using the Nikiforov-Uvarov method. 
We also obtain the normalization constants in form of the hypergeometric series. Numerical 
results show that there are only positive-energy states for bound states with spin symmetry. 
Also, the energy levels increase with the dimension and the potential range parameter a. 
Moreover, the existence of the degenerate states between E(a,n + 1,D) and E(a,n,D + 2) 
indicate that the energy levels can be completely determined using the ground state. 
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Table 1: The bound-state energy levels E n are shown in the case of spin symmetry. The 
numerical results show that the energy levels increase with both the dimensions D and the 
range parameter a. 



E(a,n,D) 








Cl — Vq — Sq - 








D 


n 


a = 0.0001 


a = 0.001 


a = 0.005 


a = 0.01 




1 


4.0032 xl0" s 


4.0326 xlO" 6 


1.0442xl0~ 4 


4.4016x10- 


-4 




2 


9.0108xl0~ 8 


9.1118xl0" 6 


2.4161xl0"" 4 


1.1130x10' 


-3 


3 


3 


1.6026xl0~ 7 


1.6271xl0" 5 


4.4842 xl0~ 4 








4 


2.5050xl0~ 7 


2.5544xl0~ 5 










5 


3.6087xl0 -7 


3.6973xl0~ 5 










1 


6.2562xl0" s 


6.3142xl0~ y 


1.6530xl0" 4 


7.1490x10" 


-4 




2 


1.2267xl0~ 7 


1.2429xl0~ 5 


3.3488 xl0~ 4 






4 


3 


2.0287xl0~ 7 


2.0641 xlO" 5 


6.0121xl0~ 4 








4 


3.0317xl0~ 7 


3.0955xl0~ 5 










5 


4.2361 xlO" 7 


4.3513xl0" 5 










1 


9.0108x10"* 


9.1118x10-° 


2.4161xl0~ 4 


1.1131x10" 


-■A 




2 


1.6026xl0~ 7 


1.6270xl0" 5 


4.4842 xlO" 4 






5 


3 


2.5050xl0" 7 


2.5541xl0" 5 










4 


3.6087xl0~ 7 


3.6973xl0" 5 










5 


4.9139xl0" 7 


5.0616xl0" 5 
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